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ABSTRACT 


Experiments have been conducted in order to determine the 
burning rates of single droplets of liquid fuels in oxidizing atmos- 
pheres. N-heptane and absolute ethyl alcohol were tested in nitrogen- 
oxygen mixtures for oxygen weight fractions between 0. 10 and 0.90. 
Measurements were made also for benzene and toluene; for oxygen 
weight fractions greater than for air, a carbon residue formed during 
the burning of these fuels, which made it impossible to obtain accurate 
data. For all nitrogen-oxygen mixtures it was found that the mass 
burning rate of single droplets of fuel is proportional to the first 
power of the drop diameter. 

The theoretical analysis on the Banaiaty of single droplets of 
fuels in an oxidizing atmosphere by Goldsmith and Penner was used to 
obtain values for the evaporation constants of the fuels tested. The 
experimentally determined values of the evaporation constants were 
found to be in good agreement with the calculated values for all fuels 


tested which did not form carbon. 
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NOMENCLATURIC 


radial distance from center of drop 

Tadius of fuel drop 

radius of combustion surface 

temperature of arnbient gas 

temperature of combustion 

weight fraction of oxygen in the oxidizing medium 


ratic of mass rate of flow of oxygen to the mass 
rate of flow of fuel vapor 


A/o Dole wo , dimensionless parameter 

thermal conductivity of fuel vapor 

binary diffusion coefficient of oxygen 

Specific heat at constant pressure of gaseous oxygen 
density of gas mixture 

density of liquid fuel 

specific latent heat of vaporization of fuel 

normal boiling point of liquid fuel 

mass rate of flow of fuel vapor 

constants in the expression (c Je =a¥ br? 
specific heat at constant pressure of the fuel vapor 
computational parameters 


evaporation constant 





I. INTRODUCTION 


During the past few years there has been increased effort to 
determine the mechanism of combustion of individual droplets of 
fuel. These basic studies were initiated in order to aid in the design 
of efficient spray injection type burners for heterogeneous combustion. 
vieterogeneous combustion is of in:portance in such widely different 
applications as stationary boilers, Diesel engines, gas turbines, and 
rocket motors. 

40st previous experi:nental work on the combustion of liquid 
droplets of fuel has been limited to the use of air as the oxidizing 


(1, 2) 


inedium, For example, G.A.E. Godsave determined experi- 
mentally the burning rates of hydro-carbon fuels such as benzene, 
ethyl alechol, and normal heptane at atmospheric sehbeece! using 
air as the oxidizing medium. Goasave estimated also the diameter of 
the flame front surrounding the liquid drop for a number of fuels. 

The mest important conclusion derived from this experimental work 
was that the mass evaporation rate of the drop is proportional to the 
first power of the drop diameter. 

A. R. Hall and J.: Died aiekoen’> have carried out experimental 
work on the combustion of drops of liquid fuels. They verified the 
result that the mass rate of burning is proportional to the first power 
of the drop diameter. Experiments were conducted also to deterrnine 
the effect of pressure on the combustion of droplets of fuel. Data were 


obtained for the burning rate at pressures up to twenty atmospheres, 


The results showed that the mass rate of burning of fuel droplets ts 
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roughly proportional to the one-fourth power of the pressure. 

in conjunction with the experimental studies, there has also been 
an effort made towards a theoretical understanding of the mechanisni 
of the cornbustion of single drons of fuel in an oxidizing ater:osphere. 


(1,2) 


Godsave obtained a useful theoretical interpretation of his ezperi- 
rental results on the assuruption that the chemical! reaction rates de 
not control the rate of burning. Uf this is the case, the proller. is 
réducéd to a study of heat and mass transfer between the flame front, 
the fuel droplet, and the surrounding atmosphere. The flame iront is 
formed at some distance from the drop surface, The oxidizer is 
delivered from the surrounding atmosphere to the region of active 
combustion by convection and diffusion; the fuel evaporates and 
diffuses, without chemical change, to the reaction front, which is 
assumed to be a spherical shell surrounding the droplet. The location 
of the reaction front i's determined empirically. Utilizing the physical 
model described above,’ lind 2) was able to derive an explicit 
expression for the burning rate of liquid fuel droplets which contained, 
however, two adjustable parameters, viz., the temperature of the 
flame front and the combustion radius. Godsave showec that the uBe 
of reasonable values ior these pararneters was consistent with the 
observed mass burning rates. 

Goldsmith and Penner (4) extended Godsave's work by using the 
hypothesis that the position of the flame front is established in such a 
way that the ratio of fuel to oxidizer corresponds to stoichiometric 
proportions. inthis way, they were able to obtain explicit relations 
for the radius of the flame front, the flame temperature, and the mass 


buraing rate, which contained no adjustable parameters. As the result 
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of an improved method of formulating the problem, Goldsmith and 
Penner were able to solve the equations without introducing invalid 
approximations concerning the physico-chemical parameters. 

There are obvious deficiencies in the assumed model for a 
burning droplet, which is shown schematically in Fig. 1. The most 
serious deficiencies are: 

1. Contrary to the assumption that a spherical reaction surface 
exists, photographs of burning droplets indicate that the luminous 
region is of the form shown in Fig. 2. Godsave measured the diameter 
of the sphere corresponding to the dotted circle in Fig. 2. The 
existing flame shape is the result of natural abe cei currents during 
burning, which have been assumed to be negligible in the theoretical 
analysis. 

2. The assumption that no chemical reactions occur in non- 
stoichiometric mixtures is obviously invalid, although a similar 
hypothesis has been used successfully in the description of gaseous 
diffusion flamed”). 

3. The calculations of mass burning rate neglect the occurrence 
of dissociation at the highest temperatures, followed by reassociation 
at lower temperatures. A more complete study, without an assumed 
reaction surface, will probably not lead to difficulties of this sort. 

Since previous experimental work on the burning of single drops 
of liquid fuels has been confined to air ao fae oxidizing medium, the 
comparison between theory and experiment is not as stringent a test 
of the theory as might be desired. With this in mind the author has 


performed experiments on a limited number of fuels in order to 
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determine the burning rate of single droplets in nitrogen-oxygen 
mixtures containing varying amounts of oxygen. The evaporatior | 
constants (burning rates) determined experimentally are compared 

in Section IV of this report with theoretical values for the evaporation 
constant obtained from the equations of Goldsmith and pees!) ° 
Although unbelievable high values are calculated for the flame tempera- 
ture, the observed burning rates are again in accord with the theoreti- 


cal results, thereby suggesting the occurrence of compensating errors 


in the analysis. 





She 
YW. THEORETICAL CALCULATIONS 
Ae Assumptions 


in order to present a clearer picture of the physical model upon 
which the theoretical analysis of Goldsmith and Pee’! is based, 
the iznportant assumptions are tabulated below. 

1. The droplets are spherical. 

2. Convection effects may be neglected. 

3. The flame front surrounding the drop is represented by a 
spherical surface concentric with the drop. All reactions ( 
take place instantaneously at this surface, at which the 
delivery rates of fuel to oxidizer are in stoichiometric 
proportions. » 

4, Steady-state conditions are assumed for fixed droplet 
sizes. This restriction greatly facilitates the mathernatical 
reatment. itis reasonable to assume that the solution 
obtained for a fixed size applies to a drop decreasing in 
size when it reaches the radius used in the steady~state 
solution. | 

5. The effect of heat transfer by radiation is neglected. 

6. otean values are used, when appropriate, for the 
physical properties. 

7. Vhe temperature of the liquid drop is assumed to be uniform 
and equal to the boiling temperature, 

8. The pressure is assumed to be uniform throughout the 


system. 
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4 schematic diagram of an evaporating and burning fuel droplet 
in an cxidizing atrnosphere is shownin Fig. 1. The radius of the 
liquid dropis Y and its temperature is the normal boiling point  . 
ahe radial distance of the cornbustion surface froin the center of the 
liquid dropletis wm and its ternperature is TO. The oxygen-inert 
gas mixture at a large distance from the combustion surface is at the 


temperature TO 


BB. Equations Necessary to Determine the Evaporation Constant 


The equations derived by Goldsmith and Penasco! ® for determin- 
ing the flame temperature T.. the combustion radius Ye , and the 


mass burning rate Mp are summarized below: 


( 
+ at -« {li+ (Yoo/ 0) ] sont i} ; (1) 
st B fi+ (Yo/ Fo)] B/% 
4 
Ye (Cro } | = _ XK « + BTo 
- = fae “ Ge dn SEE ™ 
where 


fn {14 l(Te-B/ok La + (b/2( TK + BI} 


_ap, (arbk -§)(a+ble+ §) (3) 
: (a+brh+E)(atbR-&) ” 


+e 
1 


and 
Mm, = ae. 
beef) - R/k] 
The symbols used in Equations (1) through (4) are defined in the 


Table of Nomenclature. 





me a 


The linear functional relation between my and & is 
customarily used to obtain the following expression for the variation 


} 
of droplet diameter with time: (1,2, 4) 


Peak 8% (5) 


ilere ad is the droplet diameter at time t, d is the initial droplet 
diameter, and the evavoration constant K' is defined by the following 


relation: 
R’=2 th /r % (6) 


Since K' is independent cf y , it is a convenient parameter for 
comparing burning rates of different fuels for arbitrary droplet sizes. 
The procedure for calculating the various quantities involves 
the following steps: 
(a) For suitably aabe values of the physico-chemicel 
parameters a value of T. is obtained from Eq, (1). 
(b) The computational parameter $ is next obtained frocn 
Eq. (3). 
(c) The quantity “/' is obtained from Eq. (2). 
(a) Me is determined from Eq. (4). 
(ce) Finally the evaporation constant K' is calculated from 
Eq. (6). 
Listed in Table 1 are appropriate values of the physico-chemical 


parameters used to obtain theoretical values for the evaporation 


constant. 
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Cc. Results 


The theoretical calculations for the determination of the combustion 
and temperature (T.). the radius of the reaction zone (  ), and the 
evaporation constant (K') were carried out by Goldsmith and Penner for 
combustion in air (for which the weight fraction of oxygen, Y 0,0” has 
the value 0,23). Similar calculations have been performed by the 
author for oxygen weight fractions Yo, é of 0.10, 0.37, 0.50, 0.70, 
ana 0.99, The results of these calculations are listed in Table 2, 
together with the known results for air, 

Reference to Table 2 indicates the expected trend of the data, name- 
ly, the calculated rate of burning increases with an increase in the con- 
centration of oxygen. Furthermore, with a decreased concentration of 
diluent inert gases, the combustion surface moves closer to the liquid 
drop, For a weight fraction of oxygen smaller than that for air, the 
opposite behavior is noted, i.e., for O55 = 0.10, the combustion sur- 
face moves farther from the drop and the rate of burning‘is decreased. 

The values obtained for the temperature of the combustioa zone, 
Tos are, of course, fictitious, As was stated in the introduction, the 
high flame temperatures are a direct consequence of deficiencies in the 
assumed physical model. Goldsmith and Brunet noted that proper 
allowance for dissociation and reassociation of the reaction products 
would greatly decrease the value of T. but would probably not change 
materially the temperature gradient near the droplet surface, which 
determines the burning rate. The success of the theory must be as- 


cribed to this type of compensation, at least until a satisfactory theo - 


retical model can be analyzed quantitatively. 
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Il, EXPERIWENTAL WORK 


A. Apparatus 
The apparatus used for photographing single burning drops 


suspended on a silica filament is shown in Fig. 3. The object of the 
investigations was to record the decrease in drop size with time during 
combustion. The combustion tank shown in Fig. 3 was utilized so that 
the composition of the oxidizing medium surrounding the fuel droplet 
could be varied. The tank is cylindrical, #2 inches long, with an inside 
diameter of 6-1/2 inches and an outside diameter of 8 inches. It can 
be pressurized to approximately 500 psi and can be heated to a tempera- 
ture of about 1000°F. There are a total of ten windows, 2 inches in 
diameter, that is, five windows on each side diametrically opposed. 
The tank was constructed by Mi. Goldsmith for experiments on falling 
drops, combustion avian Ftd tae, and combustion at elevated 
temperatures. 

A fuel injection system is built into the head of the combustion 
tank, This system is provided with a fuel reservoir, a fuel inlet valve 
for filling the reservoir, and with means for pressurizing the reservoir. 
The fuel system is connected through a ball check valve to a hollow 
brass needle which extends approximately 8 inches down into the tank. 

A fine silica filament was cemented to the brass needle. This 
filament was of such length that its tip extended vertically down so as to 
be in the line of sight of the top windows of the tank. It was necessary 
to thicken the end of the filament in order to retain the liquid drop on the 
filament. The diameter of the thickened end of the filarnent was approxi- 


mately 0.30 mm. 
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The drops were suspended on the filament by forcing fuel through 
the brass needle and onto the silica filafaent. The suspended drops were 
found to be between 1,5 and 1.8 mmm in diameter. As can be seen from 
Fig. 4 the suspended drops were reasonably spherical during the 
major part of their life. 

fhe suspended drops were ignited by means of an electric spark. 
A six-volt storage battery was connected through a spark coi! and an 
on-off switch to a spark plug screwed into the head of the combustion 
tank. A welding rod was used to lengthen the center pole of the spark 
plug to a length commensurate with the tip of the filament. Another rod 
was extended down from the head forming a spark gap approxiinately 
@2cminlength, across the filament tip. This method of ignition was 
found to be satisfactory for oxygen concentrations greater than about 
20 percent in oxygen-nitrogen mixtures. At Yo, e- 0,.0 it was im- 
possible to ignite a drop for the four fuels tested. 

An electrically driven Bell and iJowell "Eyemo'' 32 mm niovie 
camera was used to photograph the burning drops. The drops were 
photographed in silhouette by providing strong background illumination, 
and using a specially selected lens in order to obtain as large an image 
as possible of the dron on each frame of the film. The lens used was a 
six inch F/4. 5 Cooke Teleking Anastigmat connected to the camera 
through a 10 inch adapter tube. A lens setting of £16 was determined to 
be the most satisfactory setting. The background illumination was 
provided by a 100 watt light bulb shining directly into the tank window 
opposite that into which the camera was focused. Kodak Plus % 32 mm 


film was used. 
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The drops were photographed at a camera setting of 24 frames 
per second. The actual camera speed at this setting was checked by 
photographing a calibrated stop watch. The camera timing was found 
to be 24.60 frames per second {1 frame per 0.04065 seconds). 

A 3/32 iach ball bearing was photographed at the beginning and 
end of each 109 fvot roli of film used. This calibration was carried out 
under the same camera focusing conditions as for the burning drops, 
thereby providing an accurate reference measurement for deterinining 
the magnification of the photographs and, consequently, the actual 
size of the drops. 

The size of the burning drops was determined by using a technique 
very similar to that described by Goasave’!), The film was measured 
with the aid of a 32 mm microfilm reader (to produce iurther magnifi- 
cation) and a steel scale graduated to 1/10Gth of aninch. Two measure- 
ments were made on each frame, namely, the two perpendicular 
diameters inclined at 45° to the major and minor axes in the plane of 
observation. The mean of these two measurements was recordel as 
the "effective diameter” of the drop. If the major and minor axes do 
not differ greatly, as was the case in our tests, then it is castly shown 
that the volume of a spkere with the measured effective diameter is not 
greatly difterent tfom that of the prolate spheroid which actually 
corresponds to the shape ef the drop. 

In most cases, measurcinreats of cack drop were taken over a 
range of size extending from ignition down te approximetely une nalf 
or the original drop diameter. Tor burning in air, every fift. trame of 
the film on each Urcp was measured. Az the burning rate increased, 


with increase in oxyyen concentration, every two to four frames of the 
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film for any given drop were measured. 

In order to vary the weight fraction of the oxygen (Y 9, 3} in the 
atmosphere inside the combustion tank, a bottle of compressed oxygen 
was connected through a regulator valve to the tank. Then, for example, 
for'y > 7 8.50, the combustion tank was pressurized to a gage 
pressure of 7.1 psi. The tank was kept at the required pressure for an 
extended tirne in order to allow complete — = of the oxygen with the 
original air in the tank, The pressure in the tank was slowly decreased 
to the operating pressure, a drop of fuel suspended on the filament, the 
camera started, and the drop ignited. After each run the combustion 
tank was flushed with air before repeating the above process for another 
run. <All combustion experiments were carried out at atmospheric 
pressure. It should be noted that for Yo, + 0.10 the combustion tank 
was pressurized with nitrogen to a gage pressure of 19.9 psi before 
bleeding to atmospheric pressure. However, ignition at this reduced 
oxygen concentration was impossible and no experimental data were 


obtained, 


By. Esp erimental Results 


A aystematic study was made to determine the evaporation 
constant, <', for different oxygen concentrations and for the following 
fuels: n-heptane, absolute ethyl alcohol, benzene, and toluene, The 
experimental results gave directly the droplet diameter, d, asa 
function of the time, t. The blots of ie against t were found to be 
linear. Typical experimental plots for n-heptane and ethyl alcohol for 


Y = 0,50, and for benzene and toluene at Y = 0,23 are gsnown in 
O; © Oo, 0 
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Figs. 5to 8, respectively. Thus the linear functional relation between 
drop diameter d and time t, a’ = a’ - K't, has been verified again. 
Here d, is the diameter at time t = o (corresponding to the beginning 
of steady burning) and K' is the evaporation constant. The value of K' 
is determined directly from the slope of plots of é vs t and has the 
dimensions aa” sec '. 

The experimentally determined values of K' for the four fuels 
tested at different oxygen weight fractions are listed in Table 3. The 
results of five separate experiments were examined for each fuel at 
each value of Yo, ¥ The tabulated values of the evaporation constant 
are the mean values for each set of five runs. The reproducibility in 
each set of runs was good; the maximum spread of experimental data 
was found to be less than + ! percent. 

With the lighting used for photography during the experiments, it 
was impossible to obtain pictures of the flames suitable for measurement 
of flame surface diameter. For n-heptane and ethyl alcohol the luminous 
flame surface was not sufficiently well defined to permit conclusions 
regarding flame shape and flame size. For benzene and toluene burning 
in air, the luminous region could be seen but was not sufficiently clear 
to yield accurate measurement. However, it was possible to determine 
qualitatively that for the latter cases the apparent values of K/q@ are 
considerably smaller than the theoretical values listed in Table 2. This 
noticeable difference between the observed and calculated values for 

r/ihe has been considered in the paper by Goldsmith and Pense!”’, 


and should be ascribed to deficiencies in the theoretical treatment, which 


does not make proper allowance for the influence of convection currents 
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TABLE 3. EXPERIMENTAL VALUES OF K' (em? see” *) FOR SEVERAL 


FUELS AND WEIGHT FRACTIONS OF OXYGEN 


Weight Fractions of Oxygen 
- °° = Z 


ee 35x10 unig) ee 18, 56x 10> 


thyl 8. 58x 10 15.04%10°7 | 18.41x10> 
Alcohol 


-3 
-3 

9,85x 107? CARBON FORMATION 

7.70x107> CARBON FORMATION 


on flame shape, and which assumes a discontinuous reaction surface with- 
























out proper allowance for dissociation. 


C. Carbon Formation 

For oxygen weight fractions exceeding Yo, ae 0.23, formation of 
carbon residue was obtained during the burning of benzene and toluene. 
Photographs of toluene burning at Yo, a 0. 50 are shown in Fig. 4(b). 
Reference to Fig. 4(b) shows that the carbon formed during combustion at 
Y 6s" 0.50 remained on the filament after all of the liquid fuel had 
burned. This carbon residue was inspected and found to be of a dry, 
brittle, consistency. An effort was made to determine the value of Yo, - 
at which carbon formation is first noticeable. However, no definite con- 
clusions could be drawn since any reasonable change in the oxygen concen- 


tration above that of normal air resulted in carbon formation. The forma- 


tion of carbon in heterogeneous combustion of aromatic fuels has been 


observed by other investigators and merits additional quantitative study. 
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For oxygen weight fraction Y = 0.90, transitory formation of 


O,0 
carbon occurs during the early stages of combustion of droplets of 
benzene and toluene, but the carbon is consumed by the time combustion 
is completed, i.e., there are no deposits of the residue remaining on 
the filament at the end of burning. 

Because of carbon formation during the combustion of benzene 
and toluene, it was impossible to determine accurately the burning rates 
of these fuels. The shroudlike formation of carbon caused considerable 


distortion of the drop, and very shortly after ignition the carbon com- 


pletely surrounded and obscured the drop of fuel. 


D. Possibility of Utilizing Diameter Vs. Time Data for Studying the 





Ignition Process 


The experimental data yield not only information on steady 
burning but provide information also on the transitory ignition processes 
prior to steady burning (Fig. 9). Thus a change of burning rate takes 
aes during the time that the temperature of the liquid drop is raised 
- from its initial ambient temperature to the steady-burning temperature. 
| It is possible to study the pre-steady state burning in a manner similar 
to that employed for determining the steady burning rates a drops of 
fuel. However, spark ignition should not be used for the study of 
ignition phenomena. With the electric spark, there is considerable dis- 
tortion of the drop at the time of ignition and, in extreme cases, there 


is some loss of the liquid fuel. 
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IV. COMPARISON BETWEEN THEORY AND EXPERIMENT FOR 
THE EVAPORATION CONSTANT IN OXYGEN-ENRIC?HED 


NITROGEN-OXYGEN MIXTURES 


The calculated and observed —— constants for different 
oxygen weight fractions are contrasted in Table 4. The satisfactory 
agreement between calculated and observed values of K' suggests that 
the physical model upon which the theoretical analysis is based re- 
presents F useful first approximation for predicting the burning rates 
of single droplets of the fuels considered. The success of the theory 
must be ascribed to a nearly complete compensation between obvious 
deficiencies in the analysis: the temperature gradient near the drop 
surface which determines the burning rate, will probably be similar 
to the calculated profile if proper allowance is made for distributed 
chemical reactions, dissociation, and reassociation. 

Further comparison between theory and experiment is shown 
in Figs. 10 and 11 for n-heptane and ethyl alcohol. These plots indicate 
the calculated trend of the evaporation constant, K', (and hence of m,-) 
with increasing oxygen concentration in the oxidizing medium, and also 
the satisfactory agreement between theoretical and experimentally 


determined values. 
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V. SUGGESTIONS FOR FUTURE WORK 


Additional work on the burning of single droplets of fuel in an 
oxidizing atmosphere should include the following problems: 

1. A detailed investigation of carbon formation during 
heterogeneous burning of aromatic chemicals. 

é. coniparlact between theory and experiment for the 
effects of pressure on the burning rates of liquid fuels. 

3. Effect of temperature of the oxidizing atmosphere on 
the burning rate of liquid fuels. 

4. A detailed study of ignition phenomena. 

5. A detailed study of the effects of convection and 


turbulence on the gaseous phase of the fuel. 
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FIGURE 1, SCHEMATIC DIAGRAM OF BURNING FUEL DROP 
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FIGURE 2. SCHEMATIC DIAGRAM OF FLAME FRONT 
SURROUNDING A BURNING DROPLET 





- oe 





FIGURE 3. PHOTOGRAPH OF APPARATUS, 








Ignition Ignition t=0.936-sec. 








t= 1 71sIwere. 





tb =@9:055« 8c. 


{ 
! 
! 





SN onl" 


sec, GaOrt 5" &eer, Burning 
complete. 
4a. N-Heptane burning 4b. Toluene burning 
at Yo,o=0.50. at Yo, 07 0-90. 


FIGURE 4, PHOTOGRAPHS OF BURNING DROPS OF FUEL 
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FIGURE 5. Plot of a* Vs. t for N- Heptane and Yo, 070-50 
(K'=15.01 x 10cm’ sec.) 
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FIGURE 6. Plot of ad Vs. t for Ethyl Alcohol 
and Yojo-0.50 (K'=15.04 x 10 cm, sec!) 
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FIGURE 7. Plot of a& Vs. + for Toluene and Yo ,o= 0.23 
(K'=7.70 x 10 om.sec.) 
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FIGURE 8. Plot of a Vs. t for Benzene and Yo,o=0.23 
(K'= 9.85 x 10 cm-sec.) 
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FIGURE 9. 
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EXAMPLE OF TRANSITORY IGNITION PROCESS, 
(ETHYL ALCOHOL BURNING IN AIR) 
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FIGURE 10. COMPARISON OF THEORETICAL AND EXPERIMENTAL 
VALUES OF K' FOR N-HEPTANE. 
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FIGURE 11. COMPARISON OF THEORETICAL AND EXPERIMENTAL 
VALUES OF K'FOR ETHYL ALCOHOL. 
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